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I  -  LOW  TEMPERATURE  QUARTZ  CRYSTAL  OSCILLATOR 


GENERAL  PRESENTATION 

One  of  the  main  results  regarding  1/f  noise  obtained  during  this  program  was 
the  1/Q4  dependance  law.  A  new  senes  of  noise  measurements  was  performed  on 
different  kinds  of  resonators  to  confirm  this  law  and  principally  to  study  the 
influence  of  power  level  on  the  noise  behavior.  Comparative  measurements  were 
achieved  at  low,  medium  and  high  powers.  A  new  result  is  the  possibility  for 
guartz  resonators  to  experience  chaotic  behavior. 

Frequency  stability  at  low  temperature  was  measured  at  both  short  and  medium 
terms  using  the  dual  crystal  passive  system.  Comparisons  are  made  with  corres¬ 
ponding  stabilities  of  the  same  crystal  measured  at  room  temperature  in  the 
same  system. 

Realization  of  an  oscillator  at  low  temperature,  including  crystal  and  elec¬ 
tronics,  is  discussed,  and  preliminary  results  are  given. 


1/f  NOISE  LEVEL  VERSUS  POWER  LEVEL 
1 )  Noise  at  low  power  * 

The  measurement  system  remains  the  same  and  details  can  be  found  in  previous 
reports  (AFOSR  arants  80—01 0 S  and  81-0191). 

2.5  MHz,  5  MHz  and  10  vHz  crystals  were  measured. 


*  These  measurements  were  performed  at  the  National  Bureau  of  Standards 
F reauencv  and  Time  Division  by  0.3.  GagneDain. 


1)  5  MHz  resonators.  Six  AT-cut 
crystals  (5th  overtone)  were  stu¬ 
died.  Two  were  commercial  BVA  re¬ 
sonators  and  four  were  standard 
commercial  hiqh  quality  resona¬ 
tors.  A  typical  spectrum  (pair  of 
BVA  crystals)  is  shown  in  Fiq.  1. 

The  noise  is  essentially  1 /f  fre¬ 
quency  noise.  The  1/f3  noise  is 
due  to  the  filterinq  effect  of 
the  resonator  ;  this  occurs  at 
its  half  bandwidth.  For  three 

pairs  of  resonators  the  frequency 

1  /2 

noise  5  (f),  measured  1  Hz  from 

v 

the  carrier  were  equal  to  2.5  x 
10-13//Hz,  2.5x10“13//Hz  and 
3* 10_13//Hz.  The  raw  data  were 
diviaed  by  /2  under  the  assump¬ 
tion  that  the  noise  of  each  reso¬ 
nator  was  identical.  Obviously  by 
doinq  all  possible  pairs  one  can  obtain  a  more  accurate  value  for  each  reso¬ 
nator  independent  of  this  assumption.  All  these  resonators  bad  unloaded 
Q-factors  in  the  ranqe  2.5x106  <  Q0  <  2.7x106. 

2)  10  MHz  resonators.  Two  kinds  were  measured  :  AT-cut  and  BT-cut  crvstals. 
Both  were  third  overtone  resonators,  but  with  very  different  0-factors.  typi¬ 
cally  3xin  for  AT-cut  and  1.7xl0&  for  BT-cut.  RT-cut  crystals  have  much 
hiqhpr  0-factors  because  this  crvst al loqraphic  orientation  has  lower  internal 
losses  for  the  thickness-shear  mode.^  Noise  spectra  are  shown  in  Fiq.  1.  A 
larqe  difference  in  the  noise  levels  can  be  observed  between  the  AT  and  RT 
cuts.  At  1  Hz  from  the  carrier  these  levels  are  of  the  order  of  7x10“  3,  »Hz 
for  the  RT  resonators  to  3.2x1G“3  ,VHz  for  the  AT  resonators.  These  results 
show  a  stronq  dependence  between  the  1/f  noise  ievei  and  0-factor. 

3}  2.5  MHz  resonators.  Two  pairs  of  fifth  overtone,  AT-cut  resonators  were 

tested  with  unloaded  Q-factors  close  to  Ax  10&  .  Although  these  pairs  had 

almost  the  same  Q-factors,  the  noise  levels  differed  by  almost  one  order  of 

maqnitude.  The  presence  of  1/f3  noise  in  the  spectra  of  5V  ;f]  indicates  that 

y  o 

the  sensitivity  to  temperature  fluctuations  are  very  lame  for  this  tvpe  of 


Fig.  1 

The  solid  curve  shows  the  apparent  frequen¬ 
cy  noise  spectrum  of  2.5  MHz,  5  MHz  and  10 
MHz  resonators.  At  Fourier  frequencies  abo¬ 
ve  the  half  bandwidth,  the  spectrum  must  be 
corrected  for  the  filterinq  effect  of  the 
resonators  yielding  the  dashed  lines. 


resonator . 
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All  the  resonators  were  measured  at  powers  of  the  order  of  10  nW. 

t 

N 

The  1/f  noise  power  1  Hz  from  the  3 

w  -fit 

carrier,  evaluated  for  each  indivi¬ 
dual  crystal,  was  plotted  as  a  func¬ 
tion  of  the  unloaded  Q-factors  in 
Fia.  2.  A  linear  repression  amona 
these  experimental  points  qives 
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These  results  are  also  summarized  in 
Table  I . 


r o 
u. 


\  «  10MHZ  S’- -cut 
\  * 10MHZ  3T-cut 
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10MHz  87-out 


10 
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The  data  show  a  clear  dependence  of 
1/f  noise  on  the  resonator's  unloaded 
Q-factors,  followinq  a  1  Q*4  law.  The 
onlv  exception  is  with  the  2.5  MHz 
(#2)  crystals,  which  show  excessive 
noise  most  likelv  due  to  thermal 
transient  effects. 


io 


.'3MHZ 


.  \  9MHz 
■ 'chmz 


3VA 


\*  2 . 3MHz 

-  A 


io 


10  10 
unloaded  o-f  Atron 


f  iq. 


1/f  noise  level,  measured  at  1  Hz  from 
the  earner  as  a  funrtion  of  the  unloaded 
Q-faetor  for  the  different  resonators 
tested.  The  measurement  system  noise  is 
indicated  hv  - 


If  these  crystals  are  used  in  an  oscillator,  the  1  f  spectrum  will  qive  a 
flicker  floor  in  time  domain,  whose  cor resoondinq  values  are  qiven  in  Table 
I,  and  which  corresponds  to  the  best  stabilitv  achievable  in  that  case  with  an 
oscillator,  at  least  at  room  temperature. 


cfvstai  exit 
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2)  Noise  at  medium  power 


When  the  drive  level  of  the  crystal  is 
increased  it  exhibits  nonlinear  ef¬ 
fects  due  mainly  to  the  higher  order 
elastic  constants.  This  is  the  well 
known  amplitude  frequency  effect  where 
distortions  appear  in  the  amplitude 
and  phase  resonance  curves  and  even 
hvsteresis  can  he  developed  as  shown 
in  f lq.  3. 

The  non1 inear  behavior  of  a  resonator 
driven  1. .  transmission  can  he  repre¬ 
sented  by  the  phenomeno loqical  rela¬ 
tion 

— i  +  —  +  u2[  1  -  2e(Q)  cos  Qt] 

<1t2  Q  dt  x  i(1  +  ki2 )  =  F  cos  uit 


c  r3 
3  i 


Frequency  (Hz) 


where  i  is  the  current  through  the 
crvstal,  Q  the  loaded  Q-factor,  and  F 
the  amplitude  of  the  driving  force,  k 
is  the  nonlinear  coefficient  (related 


Fig.  3 


Amplitude  and  phase  resonance  curve 
of  quartz  resonators 


to  the  next  hiaher  order  elastic  constants'',  e  (idcosCt  introduces  a  .modula¬ 
tion  of  the  resonance  angular  frequency  w0,  which  represents  the  frequency 
noise  of  the  resonator.  Solvinq  this  equation  with  a  perturbation  method 
gives  the  phase  noise  of  the  output  signal. 


Let  Sv  (&/  be  the  f reouencv  noise  spectrum  of  the  crvstal.  When  on  vino  it  at 
low  power,  in  the  linear  range,  the  corresponding  phase  snect rum  is  given  hv 


( o ) 


Air 


At  medium  levels,  when  the  resonator  is  driven  near  the  ,;u;t  frequency,  tre 
phase  soeetrum  becomes 


Thus  the  ratio  between  the  phase 
noise  at  medium  and  low  cowers  is 

S^Chiqh  power)  (Q2  +  wo/4Q2)2 

Sjj,  ( low  power)  Q2(Q2  +  Wq/Q2) 

This  ratio  qoes  to  unity  for  Q  3>  oj0/Q 
and  is  equal  to  w2/16Q2Q2  for  Q  <K 
o)0/4Q.  Therefore  the  induced  phase 
noise  for  the  lower  Fourier  frequency 
components  can  be  qreatly  increased 
by  the  crystal  nonlineanties.  Such  a 
noise  was  experimentally  observed  on 
a  5  MHz  (5th  overtone  AT-cut)  resona¬ 
tor  driven  at  2.5  mW,  as  shown  on 
Fio.  4. 


3)  Noise  at  hiqh  power 


At  still  hiqher  power  (of  the  order 
of  a  few  watts)  the  quartz  resonator 
was  found  to  exhibit  larne  instabili¬ 
ties  as  shown  in  Fiq.  5. 


Such  phenomena  are  known  in  nonlinear 
systems  as  chaotic  behavior.  Chaos 

I 

has  been  observed  in  manv  different 
systems  ,  for  example  phase  locked 
loops2  and  Josephson  junctions 
which  can  be  considered  as  low 
0-factor  resonators.  The  main  diffe¬ 
rence  between  these  svstems  and 
auartz  crystal  resonators  stem  from 
the  hiqh  Q-factor  and  thermal 
effects,  which  stronqlv  modifv  the 
behavior  oe  the  crvstal  at  hich 


Fourier-  frequency  (Ur) 
Fig.  4 

Frequency  noise  spectrum  of  a  5  MHz 
driven  at  low  and  medium  power 


Fouriar  frequency 

F  i  a .  5 

Stable  and  unstable  states  of  a  5  MHz 
resonator  driven  at  hich  power 


s  f .  v.  -• 


Thermal  effects  will  be  predominant  mainly  if  the  crystal  is  operated  at  a 
temperature  below  its  turnover  temperature.  In  this  case  the  temperature 
coefficient  is  neqative  (of  the  order  -  4  Hz/°C).  The  temperature  rise  cue  to 
the  dissipated  power  induces  a  neqative  frequency  shift,  which  can  be  lame 
enouoh  to  pull  the  crystal  frequency  to  the  frequency  where  tne  down-jump 
phenomenon  occurs  (Fig.  3).  Thus  the  amplitude  becomes  much  smaller  decrea¬ 
sing  the  dissipated  power.  Temperature  therefore  decreases  causina  the 
frequency  to  increase  until  it  reaches 


the  second  jump  and  so  on.  This  gives 
a  cycling  with  large  amplitude  and 
nhase  perturbations,  as  shown  in  Fiq. 
6.  This  phenomenon  occurs  when  the 
amplitude-frequency  effect  and  the 
temperature  coefficient  have  opposite 
siqns. 


100ms 


Fig.  6 

Amplitude  perturbat ion  due  to  thermal 
effects  at  hiqh  power 


Above  the  turn-over  temperature  the  siqn  of  the  temperature  coefficient  be¬ 
comes  positive.  In  this  case  the  temperature  rises,  inriucinc  a  positive 
freouencv  chanqe,  which  corresponds  to  an  additional  ampliturie-frecuencv 
effect.  This  is  equivalent  to  increasing  the  nonl meant les  of  the  crvstai. 

Fiq.  !>  illustrates  the  amplitude  resonance  curve  of  a  3  MHz  (AT-cut,  fifth 
overtone)  crvstai  exriteri  with  a  7  V  rms  drivino  siqnal.  As  the  drive 
frequencv  increases  the  crystal  goes  from  stable  state  to  cnaotic  states 
through  sudden  transitions. 


This  is  quite  different  from  the  behavior  of  lew-Q  resonant  svstems.  For  tne 
latter  the  route  to  chaos  qoes  through  a  set  of  cascadinq  subharmcnic  bifur¬ 
cations.  The  rate  at  which  the  number  nf  subharmonics  increases  with  unvinq 
power  leads  quick  lv  to  a  spectrum  composed  nf  so  manv  spectral  lines  that  it 
appears  as  a  continuous  noise  spectrum. 


For  hi  ah  Q  resonators,  like  quartz  resonators,  suhharmnnir  peroration  is  net. 
observed  becausp  all  these  suhbarmonir  frequencies  are  out  of  the  resonator 
nnndwidt  h  and  therefore  filtered.  Onlv  when  their  number  is  lame  enouoh, 
such  as  when  there  are  components  in  the  resonator  harnwidth,  chaotic  beha¬ 
vior  appears.  This  explains  the  sudden  transition  to  rhaos  without  -nser na¬ 
tion  of  cascadinq  bifurcations. 


The  sideband  frequency  noise  of  the  transmitted  siqnal  was  measured  under  the 

2 

same  conditions  as  at  low  and  medium  powers.  A  f  spectrum  is  observed  for 
the  power  spectral  density  of  the  frequency  fluctuations  of  the  siqnal.  This 
corresponds  to  a  white  phase  spectrum  and  therefore  to  white  noise  for  the 
frequency  fluctuations  of  the  crystal.  This  is  in  agreement  with  other 
observations  of  chaos  ,  but  does  not  lead  to  1 /f  noise.  This  is  shown  in 
F lq.  7. 

The  level  of  this  noise  is  10-18  (power  spectral  density)  and  is  to  be  compa¬ 
red  with  the  level  of  Fiq.  1.  The  noise  is  increased  by  several  orders  of 
maqnitude. 

These  results  confirm  the  I/O4  law  of  the  1/f  noise  level,  and  this  is  shown 
quite  clearly  by  compannq  crystals  at  the  same  frequency  but  with  very 
different  Q-factors  such  as  AT  and  BT  crystals.  1/f  noise  can  be  altered  and 
maqnified  by  the  nonlinear  response  of  the  crystal  as  has  been  shown  at 
medium  power.  At  hiqh  power,  larqe  noise  is  generated,  but  this  noise  is  not 
related  with  the  previous  one,  and  the  exhibited  spectrum  is  purely  white. 


Fourier  frequency 


Fig.  7 

Frequency  fluctuations  of  a  t>  Mhz  resonator  driven  at 
illustrating  the  chaotic  state  of  the  crystal 


LOW  TEMPERATURE  DUAL  CRYSTAL  PASSIVE  SYSTEM 


The  dual  crystal  passive  system  which  is  used  allows  operation  of  a  crystal 
at  low  temperature  while  keeping  the  electronics  at  room  temperature.  This 
system  was  described  in  the  previous  report,  and  only  the  scnematic  diagram 
is  presently  recalled. 


F iq.  8 

Schematic  diagram  of  dual  crystal  passive  svstem 


This  system  was  tested  first  at  room  temperature.  When  operating  at  low  tem¬ 
perature  the  main  difference  comes  from  the  long  link,  a  coaxial  csoie, 
between  the  electronics  and  the  crystal  respectively  outside  and  inside 
dewar . 

The  phase  loop  which  locks  the  active  oscillator  to  the  passive  resonator 
acts  as  a  discriminator. 

The  overall  response  was  tested  in  real  conditions.  Quite  different  ner.avi'r 
was  observed  from  room  to  low  temperature,  as  shown  in  Fin.  da  and  °h. 

The  spurious  responses  observed  at  low  temperature  ar°  Jue  to  \fV.R  induce  ; 
mismatching  between  the  50  Q  characteristic  impedance  of  the  coaxial  earl-' 
and  the  motional  impedance  of  the  crystal  which  can  become  as  low  as  a  t«w 
ohms,  when  at  low  temperature. 


The  dynamic  thermal  behavior  of  the  crvstal  was  studied  at  mw  temperature 
The  frequencv-'-emperature  characteristic  is  represented  by  the  well  known 
relation 

Af/f  _  ,,  ' T _t  +  h  'T-T  j  +■  ...  ♦  a  dT'dt 
"  ■  ~  o  n  o  o 


wnere  ac,  bc  and  a  are  the  static  and  dynamic  temperature  coefficients. 


The  dynamic  coefficient  was 
evaluated  by  modulatinq  the 
temperature  of  the  oven  in¬ 
side  the  helium  bath  at  .1 
Hz  and  recordinq  the  corres- 
pondinq  frequency  chanqe.  As 
shown  on  Fiq.  10  the  phase 
difference  between  the  two 
siqnals  is  small.  This  indi¬ 
cates  also  a  small  dynamic 
effect.  Therefore  the  exact 
value  of  a  cannot  be  exactly 
measured,  but  its  upper  li¬ 
mit  was  evaluated  vieidinq 

a  <  6X10-10  s/K 


Fig.  iO  :  Temperature  arc  :  re  ,  ucn..  /  v  iriut ns> 


At  the  same  time  the  static 

coefficient  a  was  measured, 
o 

\ear  1.5  K  its  value  was 
a0  *  2x10-9  /K 

The  resonator  under  test  was 
a  5  MHz,  5th  overtone  AT -cut 
crystal.  Its  static  F-T  cha¬ 
racteristic  is  qiven  in 
Fiq.  11. 


Vile 


Short  term  frequency  stability  was  Treasured  bv  simultaneous  comoanson  with 
different  bic.h  stability  crystal  oscillators.  This  "tnanquiation"  method  has 
an  equivalent  resolution  which  is  of  the  order  nf  Zx’d-1"*  over  ICO  r>. 

Medium  term  frequency  drift,  i.e.  over  one  dav ,  was  obtained  bv  comparison 
with  a  Cs  clock.  It  was  not  possible  to  test  the  frequency  drift  over  more 
than  one  dav.  because  low  temperature  in  the  dewar  cannot  bp  maintained  for 
laraer  periods  of  time  without  refillinq. 

Stabilities  of  the  crystal  are  given  m  table  II  and  can  be  compared  to  the 
corresponding  stabilities  of  the  same  crystal  when  operated  at  room  tempera¬ 
ture. 


T  =  1,5  K 

T  =  300  K 

a  (10  s ) 

8  x  10-13  (BW=1KHz) 

3  x  ID"12 

drift 

V 

a  (100  s) 

3  x  TO"13 

3  x  10-12 

removed 

drift 

<  10-1 1 /day 

2x10“ 10  dav 

Table  II 

Comparative  stabilities  of  the  5  v<Fz,  5th  overtone 
AT-cut  crystal  under  test 

At  room  temperature  the  2* 10“*° /day  drift  was  observed  after  one  week  of 
oscillation.  This  value  is  not  characteristic  of  a  verv  hiph  stability  crys¬ 
tal.  But  improvement  by  more  than  one  order  of  magnitude  is  obtained  at  low 
temperature.  And  in  this  case  frequency  was  recorded  is  soon  as  the  crystal 
has  been  cooled  down,  without  any  preageinq. 

A  short  term  improvement  is  also  about  one  order  of  "scmt  ,re,  and  the  limi¬ 
tation  is  still  due  to  residual  temperature  fluctuations  acting  outside  on 


the  electronics. 


CRYOGENIC  OSCILLATOR 


Cooponents 

Different  components  were  selected  as  a  function  of  their  behavior  at  low 
temperature. 

a)  Resistors 

Metal  film  resistors  were  used.  The  low  temperature  behavior  depends  on  its 
history.  After  a  first  coolinq  (into  liquid  mfroqen,  a  drift  of  TxlfT  hou 
was  observed  with  a  30  kQ  resistor.  Temperature  cycimq  is  necessary  to  sta 
bilize  the  resistor. 


b)  Diode 

A  qermanium  diode  was  stu¬ 
died  at  1.2  °K.  The  ten¬ 
sion  Vq  which  is  0.7  V 
at  300  ° K  takes  the  value 
of  10  V  at  low  temperature 
as  shown  on  Fiq.  12. 
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Characteristics  of  r.'C  :'Ad  at 


Oscillator 


A  Colpitts-tvpe-oscillator  was  chosen  with  a  dual  qat e  MOSFET  transistor 
[7].  The  second  gate  is  used  to  control  the  excitation  level  (Fiq.  13). 
oscillator  was  operated  at  1.2  °K. 


Lid  nh  i 


Cl  =  220pF  - 


(Jjr  2  70  pF  J 


,  poiyegtrr 


|im  J:i«o  J' 


Fia.  13 


Oscillation  level  must  also  be  controlled  using  AO.' 


Fiq.  14  shows  the  dependence  of  the  output  frequency  on  the  qate  n°2  tension 
level.  In  order  to  improve  the  frequency  stabilitv,  it  is  necessary  to  reduce 
the  oscillation  level  as  much  as  possible. 


Fiq.  14 


These  results  are  preliminary.  The  oscillator  is  now  in  operation  and  fre 
quency  stabilitv  measurements  will  be  in  prorjress  in  the  near  futurp. 


GENERAL  CONCLUSIONS 


Two  aspects  were  examined  during  the  three  steps  of  this  grant 

-  the  noise  behavior  of  quartz  crystal  oscillators 

-  the  feasibility  and  performance  of  a  low  temperature  quartz  frequency 
standard. 

a)  Two  different  noise  phenomena  are  observed  in  the  fluctuations  of  the 
resonance  frequency  of  a  quartz  resonator. 

2 

-  1/f  random  walk  noise  which  is  quite  well  understood  and  can  be  attri¬ 
buted  to  temperature  fluctuations. 

-  1/f  flicker  noise,  which  is  the  main  object  of  the  study.  A  clear  cor¬ 
relation  of  1/f  noise  with  the  unloaded  Q-factor,  i.e.  the  internal  acoustic 
loss,  was  found.  The  1 /Q4  dependence  law  was  explained  theoretically  by  the 
assumption  of  fluctuations  in  the  thermal  phonon  relaxation  time  constant. 
This  important  result  can  be  a  quideline  for  the  complete  understanding  of 
1/f  noise  !  The  main  point  which  remains  to  be  demonstrated  is  a  1/f  noise 
sp;  ctrum  in  the  fluctuations  of  the  phonon  time  constant  ;  this  presupposes 
the  calculation  of  this  time  constant  which  has  not  yet  been  done. 

On  the  other  hand,  the  influence  of  impurities  was  also  studied.  It  was  shown 
that  fluctuations  of  their  relaxation  time  constants  can  also  lead  to  1/0 
dependance  law.  It  remains  to  be  shown  experimentally  that  impurity  relaxa¬ 
tion  contributes  to  the  general  1/f  noise  behavior  of  the  crystal.  This  would 
be  possible  by  measurement  of  1/f  noise  at  ^0  K  for  instance  at  the  maximum 
of  \a+  ion  relaxation. 

finally  the  noise  behavior  was  examined  at  high  powers,  when  the  crystal  ne- 
comes  stronqly  nonlinear  and  experiences  chaotic  behavior.  A  verv  iarqe  noise 
increase  was  observed.  This  noise  has  a  white  spectrum  and  apparently  ones 
not  lead  to  1/f. 


b)  The  measurements  of  frequency  stability  of  a  quartz  crvstai  resonator 
and/or  oscillator  at  low  temperature  were  mainlv  performed  by  us.inn  the  con¬ 
cept  of  a  dual  crystal  passive  svslem,  i.e.  a  passivp  quartz  resonator  m  li¬ 
quid  helium  and  a  phase  locked  oscillator  at  room  temper  at  ire.  This  supposed: 


i  h  - 


-  accurate  temperature  control  in  the  helium  bath 

-  low  noise  electronics 

-  minimization  of  cable  effects 

-  study  of  the  properties  of  quartz  at  low  temperature. 

The  results  obtained  with  a  medium  quality  quartz  crystal  exhibited  improve¬ 
ment  of  both  short  '10-IOOs)  and  medium  (1  day)  term  stabilities  by  one  order 
of  maqnitude,  and  this  without  any  preaqino. 

Stabilities  in  the  ranqe  of  a  few  parts  in  ID"  at  short  term  and  10“  per 
day  seem  to  be  feasible  with  hioh  quality  crystals.  Riohf  now  this  supposes 
principally  a  better  stability  of  the  electronics  at  room  temperature. 

The  realization  of  a  complete  oscillator  m  liauid  helium  was  undertaken.  The 
advantaqe  is  simplicity  and  the  avoid  ice  of  external  effects  at  roor  ! empe- 
rature.  The  oscillator  is  made  and  stability  measurements  are  in  pronress, 
but  results  were  not  yet  available  for  inclusion  in  this  report. 

As  a  brief  summary,  the  followinq  points  are  considered  as  the  rain  results 
of  the  program  : 

-  1 /f  noise  in  quartz  crystals  is  related  to  acoustic  attenuation  : 

-  it  follows  a  1  Q  law  ; 

-  assumption  is  made  of  fluctuations  of  thermal  phonon  relaxation  time 
constants  : 

-  impurities  can  contribute  to  1  'f  behavior  ; 

-  operatinq  a  quartz  resonator  at  low  temperature  improves  both  short  and 
medium  term  stabilities  bv  at  least  a  factor  10. 


i 
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II  -  FAST  WARM-UP  SAW  OSCILLATORS 


GENERAL  PRESENTATION 

When  a  auartz  crystal  plate  is  subject  to  a  change  of  the  ambient  temperatu¬ 
re,  it  underaoes  a  time-dependent  inhomogeneous  temperature  distribution 
which  causes  perturbation  of  the  frequency.  This  is  called  dynamical  thermal 
beha'  lor.  In  bulk  resonators,  cuts  are  found  which  are  not  sensitive  to 
varyinq  temperature. 

The  purpose  of  the  present  study  is  to  determine  the  dynamical  behavior  in 
SAW  oscillators.  In  thin  quartz  plates,  thermal  diffusion  is  deemed  predomi¬ 
nant  in  one  direction.  This  leads  to  a  one  dimensional  model  which  was  pre¬ 
sented  in  a  previous  report  (June  1982).  It  appears  however  that  in  a  finite 
plate,  diffusion  alono  another  direction  could  modify  the  tempt  ature  repre¬ 
sentation  and  the  spatiaL  repartition  of  stresses.  The  time-dependent  inhomo 
qeneous  two-dimensional  temperature  distribution  in  the  quartz  is  obtained 
from  the  uncoupled  heat  conduction  equation  subject  to  the  appropriate  boun¬ 
dary  conditions.  Since  the  heat  conduction  is  sufficiently  slow  compared  to 
the  speed  of  elastic  waves,  the  mechanical  inertia  terms  can  be  neglected  in 
the  stress  equation  of  motion,  reducing  them  to  the  quasi-static  stress  eaua 
tions  of  equilibrium.  The  time-dependent,  therma 1 lv - induced  deformation  stat 
is  obtained  from  the  static  linear  thermoelasticity. 

A  perturbation  analysis  of  the  equilibrium  equations  for  small  vibrations 
superposed  on  the  static  thermal  bias  has  been  performed. 

The  result ino  chanqes  in  wave  velocitv  show  that  ivnamic  thermal  behavior  or 
SAW  devices  is  Governed  cv  two  terms  : 

-  the  first  resulting  directlv  from  the  thermal  gradients  within  the  crvstal 
is  proportional  to  the  temperature  perturbation.  It  can  be  described  by  a 
coefficient  a  and  it  points  out  the  effect  of  a  temperatjre  uerturbat  1  on 
which  is  time-dependent  or  time-independent. 

-  the  second  one,  related  linearly  to  the  time  rate  of  change  of  measured 
temperature,  introduces  the  dynamic  thermal  coefficient  a. 


» 


THERMOELASTIC  THEORETICAL  MODELS 


1  )  Two-dimensional  temperature  distribution 


A  diagram  of  the  crystal  plate  is  shown  in  Fig.  1  along  with  the  associated 
coordinate  system. 
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Cross  section  of  the  plate 


Case  I  :  Lateral  heat  mo 
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Fig.  2 


The  xt  coordinate  axis  is  normal 
to  the  major  surfaces  of  the  pla¬ 
te.  The  length  of  the  plate  in 
the  0x3  direction  is  large  compa¬ 
red  with  the  cross  sectional  di¬ 
mensions  (2h  x  21  ' . 

We  discuss  two  special  cases  of 
heating  which  could  represent 
experimental  thermal  perturba¬ 
tions. 


Schematic  diagram  of  the  plate  heating  on 
iott'  lateral  surfaces.  Heat  transfer  occurs 
at  the  xt  z  i  h  faces 


The  first  case  was  mentioned 
in  the  previous  report.  Pres¬ 
cribed  surface  temperatures 
'  and  1'2-t'  are  at  the 
end  fares  r  *-1  arri  x;  =  -J 
IF iq.  2 ' .  A  linear  beat 
transfer  with  the  ambient  me¬ 
dium  at  zero  temperature  oc¬ 
curs  at  the  major  ‘ares  xj  = 

’  h.  Temper  it. ore  ”  ■  x  ^ .  t 
i ns i  le  t  tie  plate  is  ijoverned 
by  t  he  fol lowing  eon  at ions 
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9(xj,  x2.  0}=0  at  t  =  0  'Id' 

2 

k  =  X.2/M  is  the  ratio  of  the  thermal  conduct l v it les  in  tKe  x2  and  xi  direc¬ 
tions.  <  is  the  thermal  diffusion  constant  in  the  x2  unection  and  M  is  the 
linear  heat  transfer  coefficient. 


The  analytical  solution  of  the  diffusion  equation  (la)  along  with  the  initial 
condition  (Id)  and  boundary  conditions  (Ib-lc)  is  deduced  from  the  analytical 
solution  obtained  with  and  considered  as  time-independent,  by  Duhamel’s 
theorem  [ij.  The  result  is 


,  ,  r  cos  an  xo 

9(xi,x2,t)  =  l  2  Hh  - 1 
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with 

6  =  k2  a2  +  {in,'2%  )*  (3) 

and  an  are  solutions  of  the  transcendental  equation 

anh  tan  <znh  =  Hh  (4) 

Relation  (2)  exhibits  two  terms.  The  first  one  is  proportional  to  the  applied 
temperatures  $i't)  and  i2(t).  The  second  one  is  linearly  related  to  the  time 
rate  ,t,  and  i2(t/  cf  the  temperature  perturbation. 

Expression  (2;  has  been  calculated  for  quartz  plates  2  mm  thick  and  2  cm 
long,  a  linear  transfer  such  that  Hh  =  5Q  with  equal  applied  temperatures 
=  52(t)  =  T>  ( t ; ) . 


Figure  3  shows  the  actual  change  of  the  part  proportional  to  £>(t)  for  fixed 
values  of  X2  =  0.  and  0.9  mm  as  a  function  of  x^. 


Fig.  3 

a )  X2  =  0  mm 

b)  X2  =  0.9  mm 
Hh  =  30 

Y-cut 


rig. 


Note  that  the  temperature  decreases  very  rapidly  near  to  the  plate  ends,  so 
that  temperature  variation  is  negligible  in  the  interior  of  the  plate. 
Furthermore,  the  part  proportional  to  hit)  in  (2)  is  very  small.  Consequently 
SAW  frequency  chanqes  are  negligible  tcf.  section  II). 


A  slight  modification  of  this  problem  consists  of  maintaining  faces  at 
x 2  =  ±h  at  <4> i  ( t  )  and  £>2(t),  with  thermal  dissipation  occuring  on  lateral 
faces.  Temperature  transmitted  inside  the  plate  is.  to  a  good  approximation, 
represented  by  the  one  dimensional  model  of  an  infinite  extend  plate  with 
faces  at  ( t )  and  $2  ( t ) . 


Case  II  :  Lower  face  heatinq 
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Schematic  diaqnm  of  tne  plate  heating  nn 
the  lower  surface.  Heat  transfer  occurs  at. 
the  three  other  faces 
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The  temperature  0(xi,x2,t)  obeys  to  the  equations 
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5x2  k25x2  <  5t 
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-h  <  X2  <  h 

-a  <  xL  <  a 


{ 5a ! 


t  \i  2®  +  H6  =  Q 
5xi 


\2  —  +  H6  =  0 
5x2 


xi  =  ±  A  -h  <  x2  <  h 


x2  =  +  h  -A  <  xi  <  A 


(5b) 


(5c) 


6  =  ®(t) 


X2  =  -h  -A  <  xi  <  A 


Solving  the  equations  (5)  leads  to  the  temperature  distribution 
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an  and  u^^  are  solutions  of  the  transcendental  equations 


a  A  tan  a  A  =  HA 
n  n 
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(7) 
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(9) 


from  relation  (6),  temoeraturp  0(xi,x2,t)  can  be  written  as  the  sum  of  two 
terms,  proportional  respectively  to  *>(t)  and  b{  t)  and  such  that 

B  ( x i ,  x 2  ,  t '  =  6^  3>(t)  6^  <5(t) 


I  f  i  »i  * 
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9  and  8^  as  a  function  of  x2  are  shown  in  Figures  6.  Quantities  9^  and  9^ 
are  practically  constant  in  the  mam  ranqe  -xf  <  xj  <  x^,  x^  being  approxima¬ 
tely  8  mm  in  the  case  of  a  plate  2  cm  long. 

Moreover,  9^  is  linear  along  the  plate  thickness  and  9^  can  be  advantageously 
represented  by  a  polynomial  of  the  variable  x2  with  degree  3. 


2)  Two-dimensional  thermal  stresses  in  plates 


Thermal  stresses  arise  in  a  heated  body  because  of  the  non-uniform  tempera¬ 
ture  distribution  and  the  crystalline  anisotropy  (at  least  for  free  plates). 

A  twn-dimensinnal  thermoelastic  problem  is  considered  for  temperature  distri¬ 
butions  of  the  form  0(xi,x2,t).  For  long  plates,  the  resulting  problem  is 
plane-strain  [2] 

U1  =  ui(xi,x2)  u2  =  u2  (x^ ,  x2 )  u3  =  0  (10) 


Two-dimensional  formulations 

A  system  of  approximate  plate  equations  for  the  determination  of  thermal 
stresses  in  thin  piezoelectric  plates  is  performed  by  the  thin  plate  approxi¬ 
mation  due  to  Mindlin  [J].  Referred  to  the  0xix2X3  coordinate  system,  displa¬ 
cements  ux  are  developed  with  respect  to  x2  powers 

r  n  (n ) 

U1  =  2  x2  U1 

n=0 

v  n  ( n )  (11) 

u2  -  L  x2  u2 

n=Q 

U3  =  0 


The  n-th  order  plate  strains  take  the  following  form 


c (n)  1  '  ( n )  (n)  ,  . ,  ,,  (n+1 )  (n+1 ) ^ 

b  =  —  :  u  +  u  +  In-t-T  )  (6_.  u,  +  Ot,  u.  1 

lj  2  J-1  21  k  2k  l 


1  2  1 


and  the  linear  thermoelastic  constitutive  equations  are  written 


lj  =  CijkA  Skl  *  vij  9<xl»x2,t) 


.13. 


The  thermoelastic  constants  vXj  are  related  to  the  coefficients  of  linear 
expansion  and  the  elastic  constants  C  ^  bv  the  usual  relation 

vij  =  Ci jkJt  akJl 


The  static  form  of  Mindlin's  equations  may  be  written 


y(n) 

ij.i 

-  n 

T2j_1)  +  Fj^  =  °  n  =  0,  1,  2 

(15) 

takes 

the 

value  1  and  skips  2  and 

(n ) 

+h 

n  _ 

(16) 

'ij 

/ 

-h 

x2  T  i j  dx2 

p  (n )  _ 

n 

x_ 

+h 

■ 

(17) 

j 

2 

2  J  -h 

The  dependence  along  X3  has  been  disregarded. 


With  (12)  and  (13)  the  m-th  order  stress  resultants  take  the  form 


T(m)  =  c.  ,  ,  V  H  S<?}  -  V.  .  e(m) 
1J  ljkA  n=0  mn  kl  1J 


(18) 


H  =  2hm+n+1  /  (m+n+1) 

mn 

=  0 


m+n  even 
m+n  odd 


( 19  > 


9 )  denotes  the  quantities 


9^  ^  =  /  x 2  9(x^,X2»t)  dX2 

-h 


(20) 


The  boundary  conditions  on  the  main  surfaces 
T  2  2  ( +h )  =  0  T 1 2  ( +h )  =  0 


( 21  ) 


and  on  the  end  faces 


T'n)(*l)  =  0 
ij 


(22) 


are  that  of  a  free  stress  plate. 


• 

h 

Calculation  of  the  displacements 

a)  Basic  equations 

The  equilibrium  equations  (15)  may  be  expressed  in  terms  of  strains  bv 

«  ' 

►  ■ 

means  of  the  strain-stress  relation  (18).  The  strains  in  turn  can 
in  terms  of  displacements  u|n'  M2). 

be 

written 

i 

Doubly-rotated  cuts  have  elastic  constants  which 

introduce  coupling 

between 

• 

extension  and  flexure.  Then,  first  calculations 

are  maoe  in  the  case 

of 

sinqly-rotated  cuts,  and  elastic  constants  which 

couple  extension 

and  flexure 

m 

are  taken  to  vanish  in  the  constitutive  equations.  So,  from  (15) 

we 

may  write 

* 

the  extensional  plate  equations  in  the  form 

T'o)  =0 

11,1 

(23) 

Tvi )  _  .  o 

12,1  22 

(24) 

*  '■ 

■ 

T (2 )  _  2T  ^ 1  ^  =  0 

11,1  12 

(25) 

■ 

s 

In  these  equations,  the  boundary  conditions  (21) 

so  that  with  (17) 

have  been  taken 

into  account 

F(n)  -  0 

J 

(26) 

►.v.I; 

* 

In  the  same  manner,  the  flexural  plate  equations 

are 

•  * 

Tn.\  = 0 

(27) 

• 

T  ( 1  )  _  T  "°  )  =  0 

11,1  12 

;  28 ) 

T  ( 2 )  _  21  ^ 1  ^  =  0 

12,1  22 

(29) 

* 

The  following  notation 

« 

• 

u.,j!  3V-duj 

(30; 

has  been  used. 

*.  *,%*  J 

• 

*  ^ 

»*. 

*  * 

.v.\? 

P 

-  28  - 

t 

*.*.*_.*.  -  -  *  .  • 

V  N  . 

If  relation  '18)  and  the  constitutive  eauat ions  (12'  are  introduced 


in 


(23-251  and  '27-2q  ,  displacements  u-0',  y.^',  and  u^  on  one  hand 


are  aoverned  by  the  following  group  of  equations 


r  f  'o)  h2  (21,  n  ,M)  ,2  ( 3 )  , 
CLi  Ul,l  +  ~  ul,l  ,1  +  C12U2  +  h  U2  ' 


■  ,1 


C66  ,2h2  (3)  h2  Cl)  h14  (3),  r  , 

u,  +  —  u„  ,  +  —  u0  ,  )  1  =  C12^u 


,o. 


—  2,1  -  —  “2,1V 


hz  (2), 

—  u<  J 

3  * » * 


C22'U2 


h2  (2)l  0 

+  h  u2  1  -  - 


,o ; 


2h 


. h2  '  o )  h4  (2 )  s  _  h2  t,1)  3h  1,3; 

C,,( —  u,  ,  +  —  u  )  1  +  C  —  u  +  -  u  J  . 

H3  1,1  s  1,1,1  12  3  2,1  5  2  «  1 


e(2) 


r  h2  (2)  h2  (1  ) 

=  (V  (2  —  u,  +  —  u'  ; 
66  3  1  3  2,1 


h  ( 3 ) 
5  U2  *  1 


and  on  the  other  hand,  the  flexural  displacements  by 


r  ,  (1)  h2  h4  (3),  .  2h2  (2) 

Cll(ul,l  rUDl},1  +  Cl2  —  u2,l 


t  1  ) 

0 ,1 


2h 


=  C 


r  (  (1)  .2  (3)  (0)  h2  (2),  n 

C66  U1  +  h  U1  +  u-  ’  +  —  u-  ' -  =  0 


2,1  3  2,r 


C66  ..Co)  ^  h4  ,.(2) 


, —  u 


3  2 » * 


h  1 

—  u  +  —  u;  '  -  - 
5  2,1  3  1  5 


h  3u(?)' 
“  3ui  M 


(1) 


(3) 


2h2  (2) 


01  r  \  1  /  "  \  -/  /  \  n 

2  lc,,1 —  u,  ,  +  —  u,  ,  )  +  C  -  u 

12  3  1,1  5  1,1  22  3  2 


e;1)  1 

j 


2h 


The  matrix  notation  for  C  .  «  has  been  used. 


d)  Tree  stress  plate 


The  second  order  differential  equations  yield  terms  3vm'  as  known  expres¬ 
sions  in  xi .  Due  to  the  representation  of  the  temperature  as  a  sum  of  trigo¬ 
nometric  functions  (6),  the  solution  of  Eqs.  (31 )— (36)  under  conditions  (22. 
for  traction  free  boundaries  is  rather  complicated. 


Some  simplifications  can  be  earned  out  by  using  free  stress  conditions. 

t  \  /  o  \ 

In  order  to  allow  for  free  strain  S10,1  and  we  take  ;  3"\4l 

22  22  1  jL  J 

Z  0  and  T^}=  0  (37 

22  22 


Reported  in  (24) 


Ti"  =M 


this  condition  leads  to 


(38 


M  is  a  constant  which  vanishes  if  the  condition  ('22),  (±Jl)  =  0,  is  taKen 

into  account. 


from  (37)  and  (23),  (25)  we  obtain 


u 


(o) 

1,1 


(2) 

Ul,l 

(1) 


*  /r  * 

VC11 


22 


0(O) 

h2/5  -  e(2)/3 

8/45  h3 

e(2) 

-  e(o)  h2/3 

8/45  h5 

0(°) 

3h2/5  -  0'2) 

8/15  h3 

e  ^2 1 

-  h2 / 3 

8/15  h5 


(39 


(40 


(41 


‘.i-V 


where 


p+  P  p2  ,/p  ,  p*  —  p  ^2  yp 

11  =  11  "  1 2 7  22  ’  22  “  22  "  L12  11 


=  V1  -  V2 


C12/C22  ;  V2#  =  V2  -  V1  C12  Cli 


Flexural  equations  are  treated  in  a  similar  manner.  Then,  m  order  to  allow 
for  free  thickness  strains  S22  \  we  assume 


This  leads  to  the  flexural  plate  equations 


Tn'  =  0 


i  'll ]  =  0 

22 


T  „  p 

T 1 2  -  0 

T  (2  )  - 

T  1  2  =  0 


and  then 


'  1 )  * 
0V  ;  v 2 

- -  -y. 

4/3  hJ  (V, 


(2)  -  , 

-  U2.i  '  3 


’3  k2  J.V 

~  Jl,l 


(  1 

3  e  ■ 1 '  v* 


When  inteqration  with  respect  to  x\  is  performed  in  the  above  expressions, 
inteqrat  ion  constants  ar°  introduced  in  ° '  and  u,°  represent  inq  uniform 
translation  of  top  plate  alonq  the  xj  axis.  We  can  set  this  to  repo  without 

loss  of  qeneralitv.  Moreover,  inteuration  constants  acpenrmu  m  vp  1  and 

'  2 '  ( o )  ;  1  ' 

up  vanish  as  a  consequence  of  conditions  (21)  T  '  (t? )  =  0  and  T  „  .tl'.  =0 


As  previously  merit  toned ,  displacements,  strains  and  displacement  gradients 
are  obtained  in  a  sum  of  trigonometric  functions  which  is  not  suitable  for 
numerical  calculations. 


An  attempt  to  simplify  the  above  expressions  consists  of  writinq  the  tempera¬ 
ture  distribution  in  a  more  convenient  form  for  purposes  cf  sensitivity  cal¬ 
culations.  This  form  is  a  polynomial  obtained  by  polynomial  regression  in  two 
dimensions  (xj  and  x 2  ) . 


Polynomial  representation  of  temperature  distribution 

Charts  for  temperature  indicate  that  the  temperature  in  the  main  part  of  the 
plate  is  almost  constant  along  the  xi  direction  ;  the  region  near  the  end 
faces  can  be  replaced  by  a  polynomial  in  x^.  Consequently,  without  any  loss 
of  generality  we  may  introduce  the  somewhat  simplified  temperature  represen¬ 
tation  given  by  a  polvimn-i al ,  the  coefficients  of  which  are  given  by  a  polyno¬ 
mial  regression.  For  the  particular  case  of  interest  of  figures  (5)  (6),  9^ 
and  9*,  are  written  in  the  respective  forms 

-0.003  <  xi  <  0 

9{  =  -  500  ( X2~h '  154) 

-0.01  <  Xi  <  -0.008 

9  =  -(x  -h  )  ( -  0 . 0071  -  1 . 91  x  -  120x2';  10°  .55' 

i  2  11 

-0.0C6  <  xi  s  0 

0*  -  -  0.103  +  30.5x2  +  1.023  105x2  -  31.5  106*3  =  Fix,.  >56' 

-0.01  <.  xi  v  -0.00b 

9.5  =  Fix  '  [  5.83  +  2330  +  375.3  105x2  +  19.85  106x?[  ,57 


T he  whole  ranne  of  xi  values  is  rapidly  obtained.  M-tn  order  moments  of  9  are 
readily  performed.  Accordingly,  we  substitute  from  relations  .54  -  57'  into 
:  qs  39  -  42'  and  50  -  53  ).  Specific  calculations  have  been  made  for  guartr 
plates  2.  mm  thick  and  2.  cm  long  and  compared  to  the  icrieral  calculations 
performed  without  approximations. 


SENSITIVITY  OF  SAW  TO  TEMPERATURE  GRADIENTS 


)  Calculat ions  of  sensitivities 


When  a  quart/  resonator  •jnderqf'eq  heat  me.  two  k  mes  of  ‘.•''rerfs 
natural  frequence  :  a  direct  thermal  effect  bv  means  cf  dilatat: 
rial  temperature  coefficients  like  hq  .  ,  and  an  indirect  ther' 
feet  induced  bv  nonimeanties  of  the  quartz  crvstal. 


3  f  l~  g  r  t 
on  r?nr 

>C  C  1  ^  S  c  i 


Both  effects  can  be  considered  as  a  bias  applied  on  the  crvstal  subctra 
which  modifies  the  second  order  elastic  constants  'C  ,  beccne  ^ 

l  j  k  n  i  •  *■  - 

natural  state  coordinate  system  r  6  •  the  nonlinear  p  me  scat  ion  equation 
written 

p  u  ,  ,  =  1  A  u 
O  l,tt  lkjm  j,m  ,< 

where  p  is  the  specific  mass  and  u-  the  diso  Lace  men',  due  to  tl,e  u  i  v*  * 
quency  vibrations. 


The  boundary  conditions  corresponding  to  a  stress  rrei»  surface  are 


A  ,  u  • 

lkjm  j,m 


=  0  for 


<2=0 


A  ,  are  the  modified  elastic  constants,  which  can  be  written  m  the  f 
ik  jm 


.  .  =  C  ,  +  h  .  . 

lkjm  ik jm  lk jm 


where  the  H  ,  tensor  is  considered  as  a  small  term  with  respect  to  thi 
lkjm  K 

cond  order  elastic  constant  C  ,  It  can  be  related  to  the  temperature 

lkjm 

the  thermal  stresses  and  strains  in  the  foliowinq  wav 


H  ,  .  =  6  t,  -v  U  C,  .  +  U  .  G,  +5  C,  *  an  , 

lkjm  ij  1.0  kpmj  j.q  kima  uv  kmpquv  -lkjm 


where  t.  ,  5  and  J  are  respectively  the  thermodynamic  tensions, 

v  m 7  u .  i ,  o 

format  i  v‘s  jr  the  :i:-d lac?ment  gradients  induced  bv  temperature. 


In  the  case  c"  SAW  oscillators,  the  non  Ilf  .ear  couD.inn  cf  the  crvstal  w 
*hp  hi qn  f r**quencv  wave  depends  on  « ts  depth  pnne4,  rat  inn.  Then  perturba 
terms  h  ^ are  a  function  of  the  st°p  variable  xj  and  ran  he  written  * 
winq  the  relation 

- -  -  _  ;  c  !  n 

‘•J  -  .-V  < 

lk  im  “  "  ik ’m  ’2 


r< 

B 

% 

I 

* 


by  usinq  a  perturoation  method  5  w  6  1  ,  relative  frequency  smfts  are  caicu- 
*atea 


Au 

0 


u) 


O  O  (p)  * c  c 

r  Ad  u  i  Aq  ui 


pmP  r,k 


n!  a^;m 


n  p,q 


q  -  a_ 

p  q 


cjq  ♦  -  n 

i  —  '  a  -a  ; 
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2  p  V2  r 
o  o  “ 

p»q 


O  o 

°n. 


p. 


ui 


qP  qq 


.  0  QiS  , 

A  ,  u, 
s  k 


;  s  )  _ 

n'  and  q  are  characteristic  quantities  of  surface  acoustic  wa- 
r  s 


ves  propagating  in  a  nonperturbed  medium 


,(s) 


lu  ( t  - 

t  \  J  o 

D  V  A  0  °tSj  „ 

u,  =  /  A  u,  e 

k  L  s  k 

s 


n4  ; 


Since  the  temperature  distribution  is  proportional  partly  to  ? '  t  and  partly 

to  the  time  derivative  £(t)  (see  Eq.  6)  of  previous  section  for  example  the 

perturbation  terms  H  ,  have  two  parts 

Ik  jm  r 


1  lk  jm 


=  a 


ik  ;rn 


Jik  ;m 


Relative  frequency  shifts  Au>/<d0  are  calculated  far  i.t 


t  ions 


V  C 


and 


1 _  Au.  _  ~ 

3>(t)  wQ 


anc 


cent  nhu- 


ot- 


T ner..  when  a  SAW  oscillator  is  submitted  to  t emoerature  variations,  such  as 
fast  ‘"luetuat  ions  of  temperature  superimpose  1  on  slow  temperature  charges , 
re!  itiv->  •" require-,  sniffs  are  written 


ir^  respop  i  v  e  .  v  • 

; ? f  the  ' rrr*r.uf,r.rv  ,:a! 


tr,>?  -  r-.e 


;  f  i  '  rrerna  1  Serov  :nr  .or  d  is  f.-e  ref'  *  •••.-per  at  ,ir-*. 


Some  comments  about  coefficients  appearing  in  Eg.  67  may  be  given 

-  coefficients  aQ,  b  ,  cq  describe  the  thermal  behavior  of  the  quartz 
when  it  experiences  a  homogeneous  temperature  variation  given  by  T. 

-  a-coef f 1 ci ent  determines  the  f requencv  change  resulting  from  a  local 
temperature  perturbation  $  which  could  be  time-dependent  or  time-indepen¬ 
dent.  This  local  temperature  perturbation  induces  m  the  plate  an  mhomoqe- 
neous  temperature  repartition  followed  by  stresses  and  strains.  This  leads  to 

new  properties  of  the  deformed  crvstal  reflected  by  H  .  .  terms  (in  the  place 
r  1  k  j  m  r 

of  C  ,  .  ). 

ikjm' 

-  a-coef f lcient  is  related  to  the  time  rate  of  change  of  the  perturbation 
and  its  effect  vanishes  when  aforementioned  rates  become  neqlicnble. 


F reauencv  shifts  for  sinalv  rotated  cuts 


Coefficients  a  and  a  defined  in  relation  (66)  have  been  evaluated  in  the  fol 
lowing  cases  (see  Fig.  4) 

-  Time  varying  temperature  $(t)  is  applied  on  the  lower  surface  of  the  plate 

-  Some  linear  heat  transfer  H  occurs  on  the  other  sides  of  the  plate. 


Results  obtained  by  the  two  dimensional  model  are  shown  in  Table  1.  The  nor¬ 
malized  transfer  coefficient  JtH  is  equal  to  500. 


V,  X  AT,  X 


(ppm/K) 
(ps/K) 
ap  (ppm/K) 


a  (ps/K .) 
P 


0.096  !  0.065 


-0.13  -0.13 


ST,  X 


0.067 


-0.09 


-0.11 


Table  1 


Iheoreticnl  values  of  a  and  a  for  sinqlv  rotated  cuts 


Coefficients  and  a^  result  from  numerical  calculation  of  dynamic  thermal 
sensitivity  in  which  temperature  distribution  is  given  as  the  general  solu¬ 
tion  of  the  thermal  diffusion  equation  (see  Eq.  6). 


Coefficients  a^  and  a^  are  calculated  in  the  same  way  but  the  temperature 
distribution  is  replaced  by  the  polynomial  expansion  obtained  by  polynomial 
regression  in  two  dimensions  as  indicated  in  the  former  section  (see  Eqs. 
54-57). 


Comparison  between  a~  and  a  coefficients  on  one  hand  and  an  and  a  eoeffi- 
cients  on  the  other  hand  shows  that  for  each  cut  obtained  values  are  in  good 
agreement  and  it  seems  reasonable  to  use  a  simplified  representation  for  the 
2-dimensionai  temperature. 

It  appears  that  a-coefficients  are  weak  and  their  influence  on  the  static 
behavior  is  not  important. 


From  an  experimental  point  of  view,  measured  values  of  dynamic  coefficients 
of  SAW  devices  do  not  validate  the  adequacy  of  the  two-dimensional  model. 
Nevertheless,  we  present  a  simulation  of  frequency  temperature  curves. 


36  - 


2)  Simulation  of  the  dynamic  thermal  behavior 


From  Table  1,  frequency-temperature  characteristics  were  calculated  for 
Y-cut,  ST-cut,  AT-cut  for  some  time  rates  of  change  of  temperature. 

a)  Y-cut,  X  propagation 

— (T)  =  24.3  10“6 (T -T  )  +  18.7  10~9  (T-T  )2  +  0.091  10"6  <£(t)  -  0.11  10"6  5(t) 
f  oo 

T  =  25  °C  :  AT  =  .1  °C  and  $>(t)  =  AT  sin  Qt 
o 

®(t)  =  AT .Q  cos  Qt 

where  Q  is  related  to  the  time  varying  rate  v  of  the  temperature  by 

»  :i:2  .AT 
dt 


AT  =  .1  °C 

a)  v  =  130  K/mn 

b)  v  =  260  K/mn 

c)  v  =  390  K/mn 


Velocity  values  are 
large  and  experimen¬ 
tally  not  realis¬ 
tic.  That  means  that 
the  Y-cut  is  very 
insensitive  to  tem¬ 
perature  gradients 
compared  to  its  sta¬ 
tic  thermal  beha¬ 
vior. 


F iq.  7  :  Theoretical  dynamic  thermal  behavior  of  Y-cut. 
Oashed  line  is  the  static  frequency  temperature 
characteristic. 


b)  AT-cut,  X-propaqat ion 


Static  frequency-temperature  characteristic  of  the  AT,X-cut  presents  a 
turnover  point  at  65.2  °C.  Temperature  qradients  introduce  a  drift  of  the 
turnover  point  due  to  the  a-coef f lcient .  For  the  AT-cut,  a-coeff lcient  is 
equal  to  0.06  10-6/°K  which  would  correspond  to  a  new  turnover  point  at 
65.9  °C. 


Dynamical  thermal  behavior  of  the  AT-cut  is  computer-simulated  around  this 
new  turnover  point  for  several  temperature  variation  velocities 

—  =  3.3  10-6 (T-T  )  -  41  10-9(T-T  )2  +  0.06  10"6  *(t)  -  0.10  10“6  <§(t) 
f  o  o 

with  :  $(0  =  AT  sin  Qt 


T0  =  25  °C 
AT  =  1  °C 


a) 

V 

r  13 

K  /mn 

b) 

V 

=  26 

K/mn 

c) 

V 

=  39 

K/mn 

F iq.  8  :  Theoretical  dynamic  thermal  behavior  of  AT,X  cut 
Dotted  line  is  the  static  frequency  temperature 
characteristic  around  the  new  turnover  point  (65.9  °C) 


c)  ST-cut,  X-propaqation 


The  turnover  point  of  the  static  frequency  temperature  charactenst ic  of 
the  ST-cut  is  at  25  °C.  As  for  the  AT-cut,  the  4>(t)  term  leads  to  a  new  turn¬ 
over  point  at  a  temperature  of  25.8  °C. 

Computer  simulation  is  m^de  around  this  peculiar  temperature  value  for  seve¬ 
ral  temperature  variation  velocities  and  curves  are  plotted  in  Fiq.  9. 

Thermal  behavior  of  the  ST-cut,  X  propagation  is  described  by  the  following 
relation  takinq  into  account  both  static  and  dynamic  effects 

—  =  -39.6  1Q-9(T-T  )2  +  58.3  1 0_1 2 ( T-T  )3  +  0.067  10-6  S(t)  -  0.09  10~6  i(t- 

r  O  0 

with  :  =  AT  sin  Qt 


T0  =  25  °C 
AT  =  1  °C 


a )  v  =  1,8  K/mn 

b)  v  =  3,6  K/mn 

c)  v  =  5, A  K/mn 


F iq.  9  :  Theoretical  dynamic  thermal  behavior  of  ST,X  cut 
around  the  new  turnover  point  (25.8  °C) 


Frequency  shifts  for  a  doubly  rotated  cut  -  FST-cut 

In  the  case  of  singly -rotated  cuts,  terms  in  constitutive  equations  which 
couple  extension  and  flexure  are  vanishing.  In  doubly-rotated  cut  calcula¬ 
tions,  these  terms  are  to  be  included  in  equations.  Although  the  previous 
two-dimensional  stress  calculations  are  not  adequate  for  doubly-rotated  cuts, 
we  still  assume  that  the  assumptions  used  in  our  calculations  hold  for  the 
FST-cut  ($  =  6°  20',  0  =  -41  0  30'  and  ¥  =  26°)  and  corresponding  results  of 
sensitivities  are  presented  in  Table  2. 


FST-cut 

Calculations  by 
analytical  expressions 

(D) 

Calculations  by 
polynomial  representation 

(p) 

a  (ppm/K) 

0.17 

0.15 

a  (ps/K) 

0.06 

0.09 

To  illustrate  dynamic  thermal  behavior  of  the  FST-cut,  computer  simulation 
was  performed  around  the  new  turnover  point  a0  =  29.2  °C.  Results  are  shown 
in  Fig.  10 . 

—  =  -  20  10“9(T-T  )2  +  10  10-12 (T-T  )3  +  0.17  10-5  $(t)  +  0.06  10-6  $(t ) 
f  o  o 

with  :  <5(t )  =  AT  sin  Qt 


Fig, 10  :  Theoretical  simulation  of  dynamic  thermal  behavior 
of  the  FST-cut.  Dotted  line  is  the  static  f reguency-temoerature 
characteristic  around  the  new  turnover  point  '29.2  °C) 


In  conclusion,  the  FST-cut  exhibits  sensitivity  to  time-varying  temperatures 
of  the  same  order  of  magnitude  as  for  singly  rotated  cuts,  despite  the 
fact  that  the  model  neclects  terms  coupling  extension  and  flexure. 


GENERAL  CONCLUSION 


The  stress  state  caused  in  quartz  crystal  by  thermal  perturbation  has  been 
analyzed  for  its  effects  on  the  temperature  stability  of  SAW  devices.  A  two- 
dimensionnai  model  was  developed  in  order  to  obtain  first  ;  temperature  dis¬ 
tributions  which  are  presented  as  temperature  charts,  and  second  ;  stresses 
induced  by  thermal  effects. 

Corresponding  dynamic  sensitivity  was  obtained  by  a  perturbation  method  and 
characterized  by  two  coefficients,  the  first  one  a,  directly  related  to  ther¬ 
mal  strain  state,  the  second  one  a  reflecting  the  response  to  time  rate  chan¬ 
ges  of  the  temperature. 

Computer  simulations  of  dynamic  thermal  behavior  are  made  for  singly  rotated 
cuts  (V,  AT,  ST-cuts,  X  propagation)  and  the  FST-cut  and  for  several  values 
of  temperature  time  rates.  We  note  the  small  dynamic  thermal  sensitivity  of 
these  cuts  where  a  is  smaller  than  a  measured  for  SC-cut  bulk  wave  resona¬ 
tors.  We  also  have  discussed  the  expected  impact  of  dynamic  coefficients  on 
the  frequency-temperature  characteristics. 

It  should  be  noted  that  theoretical  a  and  a  coefficients  are  derived  in  the 
restrictive  case  of  a  free  plate.  Accordingly,  mechanical  stresses  induced  by 
the  external  mounting  have  been  ignored  in  this  treatment.  The  description 
must  be  extended  to  incorporate  effects  such  as  forces  due  to  the  differen¬ 
tial  expansion  between  mounting  and  plate.  The  specification  of  prescribed 
boundary  conditions  in  any  experimental  situation  is  a  peculiar  problem  and, 
in  many  cases,  its  solution  may  require  a  three  dimensional  model. 

At  this  point,  only  few  experimental  results  on  dynamical  temperature  effect 
in  SAW  devices  are  reported  in  the  literature.  To  confirm  the  validity  of  the 
theoretical  model,  experimental  values  for  several  cuts  .must  be  obtained. 
Consequently,  it  is  necessary  to  consider  dynamic  thermal  experimental  impli¬ 
cations  for  SAW  oscillators  built  in  the  same  technology. 

Systematic  experiments  rust  be  performed  to  move  forward  on  this  problem. 


As  a  brief  summary,  the  following  points  form  the  main  results  of  this 

program  : 

-  dynamic  thermal  behavior  of  SAW  devices  is  due  to  direct  influence  of  tem¬ 
perature  gradients  and  to  stresses  induced  by  thermoelastic  properties  in¬ 
side  the  free  plate. 

-  theoretical  dynamic  temperature  coefficients  of  SAW  devices  are  smaller 
than  those  measured  m  SC-cut  bulk  wave  resonators  (1  to  5  10~7  s/K). 

-  mechanical  stresses  due  to  the  mounting  have  to  be  included  to  model  com¬ 
pletely  the  thermal  behavior  of  SAW  oscillators. 

-  systematic  measurements  are  necessary  to  confirm  theoretical  results  obtai¬ 
ned  during  this  program. 
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